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LINE- VORTEX THEORY FOR CALCULATION OF SUPERSONIC DOWNWASH 


Bj r Harold Mirels and Rudolph C. Haefeli 


SUMMARY 

The perturbation field induced by a line vortex in a super- 
sonic stream and the doumwash behind a supersonic lifting 
surface are examined to establish approximate methods for de- 
termining the doumwash behind supersonic wings. 

Lifting-line methods are presented for calculating supersonic 
downwash. A bent lifting-line method is proposed for com- 
puting the downwash field behind swept wings. When applied 
to triangular wings with subsonic leading edges , this method 
gives results that, in general, are in good agreement with the 
exact linearized solution. An unbent lifting-line method 
( horseshoe-vortex system ) is proposed for unswept wings. This 
method is applied to determine the downwash behind rectangular 
wings with aspect ratios of 2 and 4- Excellent agreement with 
exact linearized theory is obtained for both aspect ratios by plac- 
ing the lifting line at the %-chord point. The use of lifting 
lines therefore appears promising for obtaining estimates of 
the downwash behind supersonic wings. 

INTRODUCTION 

Several methods, based on linearized theory, have been 
presented for obtaining the downwash behind supersonic 
wings. These methods utilize conical superposition (refer- 
ence 1), doublet distributions (references 2 and 3), or vortex 
distributions (references 4 and 5). Each of these methods 
has certain disadvantages. Conical superposition is re- 
stricted to wings having plan forms composed of straight-line 
segments and is cumbersome for other than trapezoidal or 
triangular plan forms. The doublet and vortex distributions 
apply to arbitrary plan forms, but provide integral expres- 
sions for downwash that are generally very tedious to eval- 
uate. The complexity of these expressions indicates that 
there is a need for a straight-forward procedure for obtaining 
reasonably accurate, if not exact, downwash solutions. 

A logical approach to the development of a simplified 
supersonic downwash theory is to derive the supersonic 
analogs of the line-vortex procedures that have proved 
valuable in subsonic theory. Certain differences exist, how- 
ever, between the properties of vortices in a supersonic 
stream and vortices in a subsonic stream. Similarly, the 
downwash fields behind subsonic and supersonic wings differ 
in certain respects. These differences must be investigated 
before an extension of subsonic techniques is possible. 


The present report, prepared at the NACA Lewis labora- 
tory during the fall of 194S, has three main objectives: 
(1) The downwash field induced by a supersonic line vortex 
of constant slope is derived and discussed; (2) the down- 
wash behind a supersonic lifting surface is examined and 
related to the downwash field induced by line vortices; and 
(3) lifting-line methods for computing downwash are pre- 
sented and calculations based on these methods are compared 
with the exact linearized solutions. Zero-thickness wings 
(lifting surfaces) are considered throughout. 

THEORY 

GENERAL VORTEX-FIELD RELATIONS 

The equations relating velocity and vorticity distributions 
in a slightly perturbed supersonic stream are derived in 
references 4 and 5. These relations are summarized in the 
following section. The velocity field is assumed to consist 
of a major supersonic free-stream velocity TJ (taken in the 
positive x-direetion) plus small perturbation velocities u, v, 
and w such that the linearized equations of motion are 
applicable. 

Three-dimensional vortex field. — The velocity field in 
vector form is 

Z=(P +u)T+ vj+wk (1) 

(All symbols used in this report are defined in appendix A.) 
The vortex vector field o> is defined as the curl of the velocity 
vector. Thus 

w=VXg 

= + + ( 2 ) 



Vortex lines are lines that are tangent at all points to the 
local vortex vector and are determined from the equation 

dx dy dz 

£ n f 
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Tlie following integral expressions for the perturbation 
velocities induced by tlie three-dimensional vortex field 
are presented in reference 4: 

- 

u ~~ 2tt f dxody 0 dz tt (4a) 

dx 0 dy 0 dz 0 (4b) 

w=-|-| f dx 0 dy 0 dz„ (4c) 

"7T I J V f o 

where 

X 0 =x— x„ 

Yo=y—yo 

z - Zq ' ■ 

/■„= iIXf—PYa*— fPZ t * 


£l r (x 0 t-zj) 

2ir Jv T 0 i 


The subscript o indicates a variable of integration. The 
integration is conducted over the volume V included in the 
forward Mach cone from the point x,y,z) that is, 


X 0 >f)jY*+Z* 


The symbol j designates the finite part of a diver- 

gent integral. (See reference 4, 6, or 7.) The procedure for 
obtaining the finite part is systematized in reference 7 as 
follows: 

The. integrals in equations (4a) to (4c) are of the form 


r_ C* M 

J^l (*2 


■A(x„) 


x 0 )»'* 


dx 0 


(5) 


The upper limit in equation (5) corresponds to limits on the 
Mach cone in equations (4a) to (4c). The finite part of this 
integral is 


IT 


= i 


’** y Yxo) 

*1 (Xi—X 0 )* i: 


dx 0 


( 6 ) 


and f are all zero, but in the z„=Q plane f=0, whereas 
£ and i? are infinite. The limits 


£'= lim £ dz 0 - 

dzQ—t 0 


7?'=lim ri dz, 

di n ->0 


are finite, however, and are given, by (reference 5) 


£'=t>B — i >T 

h'=Ut — U-b 


} 


(8) 


The perturbation velocities induced by the vortex sheet are 
ob tabled by substituting equation (8) in equations (4a) to 
(4 c). In particular, the vertical-perturbation-velocity field 
(upwash) is given by 




AV 


dx„ dy 


t> 


( 0 ) 


The area of integration -S' includes all the vorticity in the 
forward Mach cone from x,y,z. 

Line vortex. — The vortex lines through all points on an 
infinitely small closed curve bound a vortex “tube.” The 
circulation 


K=b> <r 

(where w=(£ 2 +T 7 2 +f z ) I/!1 is the resultant vorticity and <r is 
the infi n itesimal cross-sectional area of the tube) is constant 
at all points along the tube. The vortex-vector components 
at any pobit, in terms of the differential distance dl„ along the 
vortex tube, are 

dx B 

W. 


dy* 

n ^dh 

dz„ 

!=m w. 


= —J(x,)—C 

where 

0=lim[-^i-J( a:e )l (7) 

**-*iLV*» —x-o J 

The term J(x 0 ) is the indefinite bitegral of equation (5) and 
J(x i) is the indefinite integral evaluated at the lower limit. 
The justification for this procedure is presented in references 
4, 6, and 7. Finite parts are treated in a manner similar to 
that for ordinary integrals. The rules of addition, differ- 
entiation under the integral sign, transformation of vari- 
ables, and integration by parts apply. 

Vortex sheet. — If vorticity exists as a surface of velocity 
discontinuity in the z„=0 plane, outside this plane £, v, 


The elemental volume is dx 0 dy„ dz 0 —<r dl 0 . A line vortex is 
generated by allowing the cross section of the vortex tube to 
approach zero while maintaining w<r constant. The upwash 
induced by such a luxe vortex is found by substituting the 
preceding expressions into equation (4c), which yields 

w _ ff 2 1 C K(Y 0 dx e —X 9 dy,) ^ 

2tt I J c r„ 3 

where the integration is conducted along the portion of the 
line vortex within the forecone from x,y,z. A line vortex 
cannot terminate withhx a fluid flow field but must either 
form a closed curve or extend to infinity or to a boundary of 
the field. 
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UPWASH INDUCED BY LINE VORTICES 

Complicated, velocity fields can. be generated by the linear 
superposition of relatively simple fields. It will therefore 
prove useful for subsequent developments to determine the 
upwash field induced by line vortices of constant slope. 


line vortex of constant slope and strength. — The upwash 
at x,y } z due to a line vortex of constant slope m and strength k, 
intersecting the forward Mach cone at x t ,yt (fig. 1) is 
(from equation (10)) 


2 k m\y—mx) dy c 

2x Jo [(1— fi 2 m 2 )i/ e i +2m(jS i ttiy— x )y a +m i (x 2 — jS 2 y 2 — /3 2 z 2 )l 3/2 


( 11 ) 


The integration is performed in appendix B and yields 

k (y—mx^my—x) 

2ir (x i —f} 2 y i —p 1 z T ) 1/2 [(y—mx) :i +0. — f} 1 m 2 )z 2 \ 


w—- 


( 12 ) 


The finite part of the integral in equation (11) is obtained 
by substituting the lower limit into the indefinite integral 
(appendix B). No contribution appears from the upper 
limit. By a translation of coordinates, the upwash due 
to a line vortex from x u yi that intersects the forward Mach 
cone at Xb,yi (fig. 2(a)) is found to be 


w=- 


where 


k ( 7 1— mXi)(f$ 2 mT t — Xi) 

2 t rJfFj— mXi) 2 +(l— 

Xi=^x—xi 

Y 1 =y-y l 


(13) 


By superposition (fig. 2)*of a positive line vortex k from 
x u y t and a negative vortex — k from x 2 ,y- 2 , the upwash due 
to a line-vortex segment not intersecting the forward Mach 
cone is 

where the notation 

(Xi—mtXi) (£ 2 mjY t —Xt) 




rtl(T t - WiX £ ) 2 +(l— )3 2 m i 2 )Z i 2 i 
u 


(14) 

(15) 


is used. The subscript for G indicates the appropriate 
subscripts for X, 7, Z, r, and m. 

Equation (14) applies for any line vortex of constant 
strength and slope. The circulation k is positive when the' 
vortex vector is in the direction of integration. When 
the line vortex intersects the forecone, the corresponding 
limit- (infinite Gt(m) term) is neglected by application of the 
finite-part concept. 




x,x a 


M 





Figure 1.— Geometric relations for determination of upwash induced by line verier. 



(a) Line Tories from Czi.Si). 

Cb) Line vortex from fc.yi). 

(c) Line yortes from Cri,jfi) to Ca,y«3. 

Figure 2. — Superposition for obtaining upwash induced by 
line-vortex segment from (ri.yi) to On, to) . 
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Line vortex of constant slope but varying strength. — Line 
vortices may coincide to form a resultant line vortex, the cir- 
culation of which is the sum of the strengths of the superposed 
components. The general line vortex is then one of variable 
strength k along the line x 0 =x 0 (y 0 ) having the local slope 
m 0 —dy„jdx 0 . The upwash induced by the segment from - 
x x ,y x to * 2 , 2/2 may be written 


£_ KjYo-mJCo) 

2 x J fl m 0 r „ 3 


dy <> 


(16) 


infinite vortex parallel to the stream (along the lino 2/=y t ) in 
an incompressible field. 

■ These results indicate that the behavior of the infinite 
supersonic line vortex for |/Sm|>l is completely different 
from that of the incompressible-flow vortex. When jjSmj<l, 
however, both vortices have similar upwash properties in tho 
vicinity of the vortex line and are, in fact, identical for 
m= 0. _ 

2. Bent line vortices. The upwash due to a bent line 
vortex (fig. 3(c)) of constant strength is 


The evaluation of equation (16) is generally tedious. For 
the particular case of constant slope, however, this equation 
may be integrated by parts to yield a useful expression for 
upwash. Inasmuch as s 

I3 2 (Y 0 — m„X 0 ) _ d[Q 0 {m)} 
m„r / dy a 

the integration by parts gives 

w= h{\ K ° [6 ° {m)] Go{m) W dy °) (1 7) 


[Gi(m-)—G l (m+)] (20) 

where m x ~ and mfi designate the slopes of the line vortex 
before and after the bend at x x ,y x . The term r x appears as 
a factor in the denominator of equation (20) so that the 
upwash exists only in the aftercone from x u y x and is infinite 
on the cone surface (except in the 2=0 plane). This infinite 
value of upwash is not to be confused with the infinity intro- 
duced at the intersection of the line vortex with the forward 
Mach cone, which is eliminated by application of the finite- 
part. concept. 


The term K 0 [(? 0 (m)] is disregarded at a limit corresponding 
to a point on the forward Mach cone. 

Equation (17) is a generalization of equation (14) to ac- 
count for variations in strength along a line vortex of constant 
slope. Both equations are of fundamental importance be- 
cause appropriate distributions of such line vortices will be 
used to simulate a supersonic lifting surface. 

Characteristics of upwash field due to supersonic line 
vortex. — 

1. Infinite line vortices oj constant strength and slope. 
An infinite line vortex inclined supersonically (jj8mj>l) to 
the free stream is shown in figure 3(a). Because both limits 
in equation (14) are neglected, the upwash induced by this 
line vortex is zero. This result agrees with the indications 
of oblique-airfoil theory because the perturbation velocities 
are zero downstream of the envelope of the Mach cones from 
the trailing edge of a two-dimensional airfoil inclined super- 
sonically to the free stream. 

The infinite line vortex along the line (y 0 —yi)=m(x ( ,—x 1 ) 
(fig. 3(b)), inclined subsonically (|/3m|<l), has one limit that 
intersects the forward Mach cone, v r hereas the other extends 
to infinity. The upwash for 0< /3w < 1 is 


w _ L i im a (m) - « Vl-^m^- mXj ' 

2t j 2Tr(Y 1 -mX x y+(l-p 2 m z )z i 


where x lt y x is a point on the line vortex. The_ upwash is 
infinite along the line vortex. For m— 0, equation (18) 
becomes 




Y x 


(19) 


"2*- IV-j-s 2 - 

which is identical with the expression for upwash due to an 




(a) Infinite line vortex inclined supersonically to Dree stream. 

(b) Infinite line vortex inclined subsonieally to free stream. 

(c) Bent lino vortex. 

Figure 3.— Typical lino vortices. 
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LINEARIZED SUPERSONIC- WING-THEORY RELATIONS 

The perturbation velocities on supersonic lifting surfaces 
(zero-tliickness wings), have been evaluated for a large 
variety of plan forms. (See, for example, references 7 to 9.) 
Those results will now be utilized to determine the vortex 
field generated by a lifting surface. 

Velocity potential. — If the boundary conditions for a lifting 
surface are specified in the 2=0 plane, the u and v velocities 
are antisymmetric and the w velocities are symmetric about 
this plane. Thus the velocities on the top and bottom sur- 
faces of the 2=0 plane are related by 



u T = —v . B ' 
l' t = — d b >- 
W T =W B , 


( 21 ) 


The discontinuities in u and v constitute a vortex sheet. 
Because the Sow is everywhere irrotational, except across 
this sheet, a perturbation-velocity potential 4> can be so de- 
fined .that 


d<j>=^dx+^dy±^dz ^ 
=u dx + v dy+w dz 


( 22 ) 


The undisturbed flow field ahead of the wing is considered 
to be of zero potential. The boundary conditions require 
that u T —u B — 0 off the wing (antisymmetry of u and zero 
lift off the wing). The potential in the 2=0 plane may then 
be obtained by integrating along lines of constant y 


<t>T — 

( u T dx 

J 

lx t 


(\E 

4>b— 

v 

u B dx 

<Xi 


(23) 


where x t is the equation of the leading edge as a function of y. 
From equation (23) it may be concluded that 

1. Everywhere in the 2=0 plane except behind the wing 
leading edge <j> r =<p B —0. 

2. At a given span station, <f> T and <j> B remain constant for 
all values of x downstream of the trailing edge. 

Lines of constant potential for rectangular and triangular 
wings are shown in figure 4. 

Vortex lines. — The equation for the vortex lines (from 
equation (3)) is 


i)' dx—£' dy = 0 


(24) 


When the following expressions (from equations (8), (21), 
and (22)) are substituted 


% r — — (p t — »b)= — 2 


by by 


■t}' = u t —u b =2 


r g dfe 

br br 



(a) Rectangular wing. 

(b) Triangular wing. 


Figure 4— Lines of constant potential for rectangular and triangular wings. 

the equation for the vortex lines becomes 

tt *+! py-w dl+ w dv=0 C25> 

Comparison with equation (22) shows that equation (25) 
represents lines of constant potential. Thus the vortex 
lines coincide with the lines of constant potential in the 2=0 
plane. 

Circulation. — The circulation included between two points 
x u yi and Xi,yi on a wing is given by 

K== tj) u dx + D dy-\-wdz (26) 

The path of integration is arbitrary except that the path 
should cross the 2=0 plane only at the two specified points. 
If the integral is taken along the top and bottom surfaces 
of the 2=0 plane, 

i 

C afe, *2 f*lrVl 

{u T dx+v T dy)+ ( u B dx + c B dy ) 

l,fl Jz 2.r2 

= (®r, 2 — 4> t, 1 ) + (4>b, i — ©b, 2 ) = 2 — A<£i (27) 

where A<£ equals 4> t —<I>b and represents the jump in potential 
at the point. The quantity A<f> is, in fact, the doublet 
strength (reference 2), so that the net circulation between 
two points equals the difference in the doublet strength 
between those points. The equivalence of a doublet dis- 
tribution and a vortex distribution indicates that the flow 
about a lifting surface can be calculated on either basis. 


9*j 664(5 — ol- 


42 
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Circulation and lift. — The lift per unit span is given by the 
chordwise integration 


consist of elemental vortices of infinite length along 
2 /=constant lines, each having the strength 


L{y)=\jU*^ (C P ,B-C,. T ) dx 


(28) 


dK=? dy 0 = — 


dT 

dy 0 


dy 0 . The vertical perturbation velocity at 


After substitution of the linearized values for pressure co- 
efficient 


C , 


P,T- 


2u r 

U' 


2 0f/>7' 

VHx 


os ,y,z due to the elemental vortex along y—y„ is, from equa- 
tion (19), 



n — 2 u b_ 2 d<t> B 

Lp - B TT~ TT~Zx 

and integration, equation (28) becomes, because A0;=O, 
L{y)—pU A 4>t (29) 


so that the vertical-perturbation-velocity field at infinity is 
given by 


10 = 


1 fi Y, dT ,, 

2 xj_5 Y?+z* dy, y ° 

2 


(33) 


The factor A 4> t is the circulation included between the leading 
and trailing edges at the spamvise station under consideration. 
When this circulation is designated r, equation (29) becomes 

L(y)= P Ur .. (30) 

which is the familiar incompressible-flow relation. Also 

dr d(<j>T — 4>b)i _ 

dy dy 


where 6/2 is the semispan. The velocity field obtained from 
equation (33) is identical to that induced by a subsonic wing 
with the given span loading. The velocity field at infinity 
is thus independent of Mach number (excluding the influence 
of Mach number on dT/dy a ). This result has been derived in 
references 1 and 2 by other methods. The evaluation of 
equation (33) is relatively simple and may be used to approxi- 
mate the downwash several chords behind a supersonic wing. 

REGIONS OF INFINITE DOWNWASH 


= (flr— v B )t 

(3D 

Equation (31) relates the shed vorticity to the rate of change 
of spamvise loading. 

APPLICATIONS TO CALCULATIONS OF SUPERSONIC 
DOWNWASH 

The vertical pertubation velocities due to a supersonic line 
vortex of constant slope have been presented in a form that 
permits analytical or mechanical evaluation (equations (14) 
and (17)). The vortex distribution associated with a 
lifting surface has also been discussed. These relations will 
be used to develop exact and approximate methods for calcu- 
lating downwash behind lifting surfaces at supersonic speeds. 

DOWNWASH AN INFINITE DISTANCE BEHIND WING 

The vertical-perturbation-velocity field behind a super- 
sonic lifting surface (from equation (9)) is 


w= — 


g‘ir (Y^'-xYn 
2 v J r„ 3 


d x„ 



(32) 


where the integration over the plan form is designated by 
S p and over the wake by S w . As x becomes infinite, X 0 
also becomes infinite in the integral for the bound (plan- 
form) vortices. This integral then becomes zero because 
X 0 is of higher order in the denominator than in the numer- 
ator. Thus, only the integration over the trailing vortex 
sheet contributes to the vertical perturbation velocities at 
infinity. The trailing vortex sheet may be considered to 


Approximate solutions may modify or introduce singu- 
larities in the downwash field. It is therefore of interest to 
establish the regions for which linearized theory indicates an 
infinite downwash. 

Infinite downwash in z=0 plane.— The vertical perturba- 
tion velocities in the 2=0 plane that exist an infinite distance 
behind a wing having a discontinuity hi dT/dy , will first be 
considered. These velocities are determined from the 
equation 


b 



The Cauchy principal value is required for points on the 
vortex sheet. For the particular case of triangular loading 
(fig. 5(a)), 

dr 2r„ 

dy* b 

in the integration interval— ^<y„<0 and 


dT 2T m 
dy„ . 6 


in the interval 0 <y t < ^ (F w is the circulation at the miclspan). 


When these values for dTjdy 0 are substituted in equation (34), 
the integral yields- 




r— 


( 35 ) 
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(a) Triangular load distribution. 

(b) Elliptic load distribution. 

Figcbe 5.— Upwash in z=0 plane an infinite distance behind 
Wings o£ triangular and elliptic loading. 

Infinite upwash exists along the lines y— ±b/2; whereas 
infinite downwash exists along the fine y= 0. These infinite 
values are due to the discontinuities in the span vise vorticity 
distribution and apply for all points on these lines down- 
stream of the wing trailing edge. 

In general, if (dT/dy a )~ and {dT[dy 0 ) + represent a dis- 
continuity in the rate of change of spanwise loading at 
station y=yu infinite vertical perturbation velocities will 
exist along y~yi downstream of the trailing edge. For 
(drjdy a )~<^(dr[dy 0 ) + infinite upwash will exist, and for 
(dTjdy„)~^>(drldy e ) + infinite downwash will exist along this 
line. Such discontinuities in dF/dy e originate both at a 
wing tip and at those points along a supersonic trailing 
edge where the plan-form slope is discontinuous and the 
local-wing-pressure coefficient is not equal to zero. This 
discontinuity in vorticity may be verified by the methods 
employed in reference 1 for finding the upwash and the side- 
wash directly behind a supersonic trailing edge. 

The discontinuity in shed vorticity at the tips of an ellip- 
tically loaded wing is a special case of the previously stated 
rule governing infinite vertical perturbation velocities in 
the z=0 plane. For wing loading given by 

r=r ^/i— pV? 


Integration yields 


for 
and 

w 
for 

|y|>! 

As indicated in figure 5(b), the vertical perturbation ve- 
locity is discontinuous at y— ±6/2, but is bounded for all 
points on the vortex sheet. 

Infinite downwash on Mach cone from wing tip. — Refer- 
ence 3 indicates infinite vertical perturbation velocities on 
the downstream Mach cones from the tips of a triangular 
wing. This result will now be extended to apply to any 
wing tip formed by the intersection of a subsonic leading 
edge and a supersonic trailing edge, provided that the slope 
of the subsonic edge is not zero at the tip (fig. 6). 

The contribution to the vertical-perturbation-velocity 
field due to the bending of an elemental vortex at the trailing 
edge (fig. 6) is, from equation (20), 

dr , 

w dy ° 

dw= [<?. (m -)—<?„ (0)] (36) 

Air 



where 


O,(0)=- 


•1 X.Y, 


r a (F/4-2, 2 ) 


The vertical perturbation velocity at a point on the Mach 
cone from this tip, due to the bending of the vortices, is 
found by integrating (along the trailing edge) 




(37) 
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r 



Figure 6.— Wing tip formed by intersection of subsonic leading 
and supersonic trailing plan-form edges. 


Equation (37) in the expanded form becomes 


A lifting line concentrates the ehonhvise loading into a line. 
Thus the bound circulation is represented by a line vortex 
of variable strength (K=r=A$L,); whereas the trailing- 
vortex. system maintains the same vortieity (£' —~dTjdy„) 
but now originates at the line rather than at the trailing edge. 

Bent lifting line. — A lifting line approximating the section 
centers of pressure 

f * (C v , B —C r >, T )(x— r,) dx 

x cn=Xri 

(C'p.b — C r p,r) dz 

seems to be a reasonable representation of a swept bark or 
sweptforward wing. In order to facilitate dovvmvasli 
calculations, the line of section centers of pressure can be 
approximated by two straight-line segments, each connecting 
an end point to the midpoint of the lino. The result is a 
bent lifting line (fig. 7) of span b and effeetivo chord c. 
(r-distance between midpoint and end points of lifting line). 

I V 


J_ f° ? [Y-imnXJW.MYo-Xj j 
2 TrJ fa \[Y 0 -(m 0 -) A r fl ] 2 +{l-j3 2 (mr)^ 2 

XqYq l 1 dT , 

IV+s 2 } r, dy 0 ay ° 


(38) 


The limits of integration are roots of r„, so that the factor ljr 0 
is singular of order }i at the limits and the integral is 
improper (assuming 0 ^0). The convergence, of 

the integral depends on the nature of dT/dy 0 at the limits. 
Appendix C shows, however, that dTjdy 0 is also singular of 
order % at the tip. The combined singularity causes 
equation (38) to diverge at the upper limit and results in 
logarithmically infinite vertical perturbation velocities on 
the. Mach cone from the tip. The divergence, is a conse- 
quence of both the singularity in dVjdy 0 and the singularity on 
the Mach cone from a bend in the elemental line vortex. 
The infinite vertical perturbation velocities on the Mach 
cone from the tip do not appear in the 2=0 plane inasmuch 
as equation (38) then reduces to 


w— 


2 tr 



V Xp~ $ 2 YJ dT j 
[Y-{mnX 0 ] Y 0 dy 0 ay ° 



Tiguke 7.— Bent-Hfting-lIoe approximation lor swnithuck wing. 


For a sweptback lifting line, the slopes at the midpoint are 


^ 2c, 

and 

( m o + )t a ~o~2^ 


and the singularity due to r 0 is no longer present. 

APPROXIMATE DOWN WASH SOLUTIONS 

Several approximate methods for obtaining downwash 
were considered. Methods based on a lifting line seem the 
most promising and are discussed in the following sections. 
A bent lifting line is proposed for determining the downwash 
behind a swept wing and an unbent lifting line (horseshoe- 
vortex system) for determining the downwash behind an 
unswept wing. These methods are applied to compute the 
downwash behind triangular and rectangular wings and 
the results are compared with the exact linearized solutions. 


The vertical perturbation velocities induced by the bound 
vortices can be determined from equation (17) and yield 

— G'{m+)] V ' G 0 {m) ~ dy 0 (39) 

where the integration is conducted along the lifting line 
(figs. 8(a) and 8(b)). The vertical perturbation velocities 
induced by the trailing-vortex system are obtained from 

GM w. iy ‘ (40) 
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The term <?„(0) is as defined for equation (36) and the 
integration is again conducted along the lifting line. The 
vertical perturbation velocities induced by the complete 
lifting-line system (bound and trailing vortices) are then 
the sum of equations (39) and (40) and yield 

»«£=■ [G 0 (m~)- <?,(m+) P 4 mrn)-G o m dy 0 

■l-k ° atj ay 0 

(41) 

The value m 0 =(m,,~) ¥ =0 is used for the integration interval 
V<l<Vo < 0 and m 0 ={m,f) tlo= a is used for the interval 
0 <y a <y<,. 

Equation (41) is the expression for the vertical perturba- 
tion velocities behind a swept wing using the bent-lifting-line 
approximation. This equation can also be considered as 
derived from a superposition of a system of bent line vortices 
of constant strength (equation (20)), as indicated in figure 9. 
The nonintegral term of equation (41) is then the contribu- 
tion of the bent line vortex of strength whereas the 
integral term represents the contribution of the bent elemental 

line vortices of strength 


, dr , 

± dy 0 dy ‘“ 


The integral term of equation (41) may be evaluated by 
analytical or mechanical methods. When mechanical 
methods are used, the singularities in the integrand must 
be isolated. Suitable procedures for isolating the commonly 
encountered singularities are as follows: 

(1) Singularity due to intersection of forward Mach cone 
with lifting line. The integrand in equation (41) is infinite 
at the intersection of the forward Mach cone and the lifting 
line (for 2^0). The contribution Sw of the integral, for the 
interval 2/4— (fig. 8(a)), to the vertical-perturbation- 
velocity field may be written 




l( c '’ a) 


(aj Intersection of forecone with lifting line. 

(b) Intersection of forecone with edge of vortex sheet, 
(c) Point on vortex sheet. 

Figetre 8.— Improper intervals of downwash integral. 


where e* is a convenient length. The first integral in equation (42) is proper and may be mechanically evaluated. The 
second integral is, for 0<jSm„<(l, 

f ** dy 0 ( m Q . r , ,)n 

To ~ [ m «Vl-^ 2 m 0 2 r 0 -{-(l-/3 2 ffl a 2 )2/«-m 0 (r-/3 2 m 0 2/)jj^_^ 




Vl -Pm* 


sinh 


/ c 

V 2[m 9 (r — 


(1— g 2 m a 2 ) t b 




for /S m,> 1 , 


P r ‘ dy °J\ m ° ! m o( x — i6 2 m,y)— (1— j5 2 m 0 2 )i/,, "j i ' i 2m, . / 

r o l_y/3 2 m 8 2 — 1 A 8m 0 \(y—m 0 x) 2 (l—p 2 m<r)z 2 Vj3 2 m„ 2 — 1 V 2[m 0 (x— 

f '* dyo r m 0 \2(y—m 0 x)y a 

Jt s-«6 r B L y— i 


Q3 2 mf~ l)t a 


|3*m # y)— (1— 


and for |3m 0 = 1, 


a +m 0 3 x i —y 2 —z i 
■ rn a x 


n OmL 

n-«4 m °\ m„ x— 


where m 0 is the slope of the lifting line at y b . A similar pro- 
cedure applies for the singularity at y a . The contribution 
Sw of the interval y^y a <?/«+«« to the downwash at x,y,z 
may be computed from the preceding equations by reflecting 
the lifting line about the .e-axis and then computing the con- 


tribution of that interval to the downwash at x,—y,z. 

(2) Singularity in dYjdy 0 at wing tip. The vertical per- 
turbation velocity at a point the forward Mach coBe of which 
intersects the edges of the tra iling vortex sheet (fig. 8(b)) is 
obtained from equation (41) by integration between the 
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Figuke 9.— Superposition of bent line vortices yielding bent-lifting-line representation of 
suporsonic wing (equation (41)). 


limits <y,< 7 p 


If dT/dy a is singular at these limits, the 


singularity may be isolated by a procedure similar to that 
used in equation (42). Thus, the contribution to the vertical - 
perturbation-velocity field, from the integration interval 

6b<y 0 < may be written 


^ dy 0 

2 ~ tb < 

= ~^/ 6 2 ([<?o(»l)-^(0)]- 

T _e * 

{[G 0 {m)-G o m} v ^^ o dyo- 

Jl 

T, | \dy. 

2 6 * 

where 


(43) 


x ; 


I dv, _ 

j„. dy 0 - — (r)_, _6 


— «6 


dy„ 




A similar procedure applies for the singularity aty„=— ~ 

(3) Singularity at points on vortex sheet. The Cauchy 
principal value of the integral (equation (41)) is required 
for determining the vertical perturbation velocity at a point 

on the vortex sheet z = Q,~ \<y<^- For this case, when 

the interval y— e<y 0 <y-\-e (fig. 8(c)) is considered, the 


contribution of the integral to the vertical-perturbation- 
velocity field is 


Sw 


1 f»+« m^ xf—pY} dT 
m 0 X 0 ) Y„ dy. 


* I ° v 

(1« — 


dy. 


(44) 


If dTjdy 0 can be approximated by the first two terms of the 
Taylor’s expansion 


dr 
dy o 


v 0 ~v 


+ (Vc 




dyffi.-v 


and e is sufficiently small that 


mAXf-pY'A. 

( Y,-m,X, ) “ 


equation (44) may be written 




(45) 


Inasmuch as dTjdy, is an odd function in y, for a wing 
symmetrical about the midspan, dTjdy, is either discon- 
tinuous or zero for y„—0. If dr/dy 0 is discontinuous, the 
vertical perturbation velocities are infinite along this line. 

Unbent lifting line. — The unbent lifting lino (horseshoe- 
vortex system) appears to be a reasonable representation for 
an unswept wing. The use of such a lifting line considerably 
simplifies the calculations. 

The vertical perturbation velocities induced by an unbent 
lifting line (along the y-axis) are obtained from equation (41) 
by setting m,— » and are equal to 


1 C Vb xYoW-pz*) dT , 
2 tt Jy 0 r 0 (x J -M(IV+ 2 *) dy. ay ° 


(46) 


The singularities are isolated as follows: 

(1) Singularity due to intersection of forward Mach cone 
with lifting line. 


8 w— 


_ir r * r zY.w-tJ’z*) dr 

2ir Jn-H L(ri-,8 i A( r o i +^) dy, 

7 £VdT\ "1 dy, 1_ / fPzf _dT \ f»‘ 

\xY 0 dy,)y,=y J r. 2x \ r Y, dyJi t m,J n ~ ti 


dy : 

r. 


where 


X 


'* &_|sin 

Vi— e t r o ft 



(2) Singularity in dT/dy , at wing tip. 


J_f 2 ( xY.{r.*-pe*) r ” 1 ) dr 

2 rji {r,(x*-/3*z*>(Y.*+z 2 ) Lr,(z 1 -j3 2 z 2 )(Y/+z 2 )l b \Wo 
T _e ‘ ( '•tJ 


, , 1 f xYlrf-Pz 2 ) 1 

ay ° ir 2TLr <) (x*-!5 2 z%YS-hz*)J 4 

Vo"“ 2 


(F> JL_«, 

Vo *» ** 


(3) Singularity at points on vortex sheet. 



The chordwise location of the unbent lifting line that will give 
the best average agreement with the exact linearized solution 
is still to be determined. According to the techniques used 
in subsonic-wing theory, an unbent lifting line at tlio wing 
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center of pressure should be a good first approximation; 
however, further investigation is required. It may be pos- 
sible to determine the best location for each general class of 
plan forms by comparing the lifting-line solution with the 
exact linearized solution for several representative plan forms. 

Examples and discussion. — A bent lifting line and an un- 
bent lifting line will be used to obtain solutions for the down- 
wash in the y — 0 plane behind triangular wings having sub- 
sonic leading edges. The ehordwise distribution of wing 
loading suggests the use of the bent lifting line. The unbent 
lifting line will be used for purposes of comparison. An 
unbent lifting line also will be used to compute the downwash 
in the y=Q plane behind rectangular wings. The ehordwise 
location of the unbent lifting line giving the best average 
agreement with the exact linearized solutions will be noted 
for both the triangular and the rectangular wings. In order 
to simplify the expressions, Af=y2 (that, is, (3=1) will be 


assumed. 

The spanwise circulation distribution for a triangular 


wing of span b and root chord c t , having subsonic 



leading edges, is (from equation (24) of reference 2) 


r=A^ 


aUb r 

= ~E~\ 



where E is the complete elliptic integral of the second kind 


with modulus 
is 


V i- (4 


)• 


The circulation at the midspan 


T m = 


aUb 

E 


and the rate of change of circulation is 

dr _ 4 aU y„ 

d Vo~ Eb / 4y.* 


The expression for the vertical perturbation velocities in the 
y = 0 plane, using the bent-lifting-line approximation of 
figure 10(a), is 


aUb rri , , 4cd7 f 

W ~ tE M r.-o+ T 6^J 0 


iaU [g.(m+)~ Q.(Q)1 


V 1 b iVo 


Vo dy 0 (47) 


where (m, + ) ra=0 =~ Equation (47) is obtained from equa- 


tion (41) and the properties of an even function. The vertical- 
perturbation-velocity field b ehin d an unbent lifting line 
having the same loading (fig. 10(b)) is, from equation (46), 


4aj7 rn rOy 8 — g *)y* dy 0 

TbEja I Z 

r„(x 2 - g 2 ) (y *+ Oy 1 -p Vo 1 


( 48 ) 






y,y° 


A**t 


V 


*«. ..--Unbent lifting line 
-- \ 

WK 




(c) 

(a) Bent-lifting-line representation of triangular wing. 

(b) Unbent-lifting-line representation of triangular wing. 

Cc) Unbent-lifting-line representation of rectangular wing. 

Figure 10 -Lifting-line representations of triangular and rectangular wings. 

i 

The spanwise circulation for a rectangular wing of aspect 
ratio — > 2 for the tip region \—c T <y„< xis 

C T Z 6 


r=A<£ ( = 


4aU 


hj (yo + e,— f ) (| -y 0 )+ Cr tan 1 I 


n 


2~ V * 


and for 0<y„< c T 


r= 2 aU c T 


Vo+Cr— J 
(49a) 

(49b) 


Equations (49a) and (49b) were obtained from equation (20) 
of reference 8. The wing is illustrated in figure 10(c). The 

rate of change of spanwise circulation for ~—c r <y a < ^ is 


dr 4 aUc, 

dy — 


I 


b 

VoT^r „ 


' M I- 


Vo 
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The vertical-perturbation-velocity field in the y=0 plane 
behind an unbent lifting line having this loading is 


x(.r»-s*)y a '~\ 

^-zW+z*)- 


/ , b 

J yo~\~Cr 2 


Equation (47) was evaluated to determine the downwash 
along the line y= 0, z — 0 and along the line y=0, z= 0.1 b 
for triangular wings with aspect ratios (2 bfc T ) of 1.6 and 3.2. 
Equation (48) was. evaluated for the same wings at points 
along the line y=0, z— 0. Equation (50) was evaluated at 
points along the line y= 0 , 2=0 for rectangular wings with 
aspect ratios (b/c r ) of 2 and 4. The integrations were 
mechanically performed for equations (47) and (50) and 



1.0 1.2 /.4 1.6 /.a. _ 2.0 2.2 

Chord behind apex, — 


(a) t=0 plane. 

(b) i=0.I b plane. 

Fiocee XI.— Downwash in i/=0 plane behind triangular wing with aspect ratio 26/Cr of 1.6 
using bent lifting line. M=-^2. 



Chords behind opex, 


Ca) 1=0 plane. 

Cb) *= 0 . 1 6 piano. 

FXQOSE 12.— Downwash in i/=0 piano behind triangular wing with aspect ratio 2 b,'c, of 3.2 
using bent lifting line. 

analytically for equation (48). The results are compared 
with the exact linearized solutions obtained from references 1 
to 3 in figures 11 to 14. 

The bent-lifting-line solutions for the triangular wings are 
shown in figures 11 and 12. The discontinuity in the curves 
for the 2=0.1 b plane (figs. 11(b) and 12(b)) designates 
large negative (upwash) values that become infinite on the 
Mach cone from the tip, as indicated in the discussion asso- 
ciated with equation (38). The agreement with the exact 
solutions is good for all points except those within % chord 
from the trailing edge. The correlation is unexpected 
because .of the large contribution of the noninlegral term of 
equation (47). This term is associated with the bending of 
the lifting line at the midpoint. This bend was artificially 
introduced. The agreement with the exact downwash 
solutions, however, indicates that the bent lifting line is a good 
average representation of the triangular-type wing (at least 
in regard to the downwash solution in the region of the line 
y= 0, 2=Q). This method should give even better agree- 
ment when used to represent a sweptback wing such as the 
one represented in figure 7, because the bent lifting line would 
then more closely approximate the actual vorticity distribu- 
tion. 

The unbent-lifting-line solution for the triangular wings is 
presented in figure 13. The lifting line is placed at the 
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? 4 -cliord point to give the best average agreement with the 
exact linearized solution. (The center of pressure for these 
wings is at the % -chord point.) The agreement is not as 
good as that obtained with the bent lifting line and indicates 
that the bent line is more suitable for computing downwash 
behind triangular wings. The smaller-aspeet-ratio wing 
(26/c r = 1.6) is in better agreement with the exact linearized 
solution because of the rapidity with which the downwash 
approaches the asymptotic value at- infinity. 

The unbent-lifting-line solution for the rectangular wings 
is presented in figure 14. The downwash obtained with the 
lifting line at the J4-chord point is in excellent agreement 
with the exact linearized solution. (The location of the 
center of pressure is at the /4-chord and %-chord points for 
aspect ratios of 2 and 4, respectively.) These calculations 
suggest that the unbent lifting line is well suited for com- 
puting the downwash behind unswept supersonic wings and 



Figure 13.— Downwash along line y=0 r z=0 behind triangular wing using unbent (horse- 
shoe) lifting line at ®£-chord point. M=-yt «T 



Figure 14. — Downwash along line g=0 r z=0 behind rectangular wing using unbent (horse- 
shoe) lilting line at JS-chord point. 2. 


' that the best ehordwise position for the unbent lifting line 
is at or slightly downstream of the wing center of pressure. 

A more accurate estimate of the downwash in the z = 0 
plane close to the trailing edge maybe obtained by judiciously 
fairing the curve obtained by the line-vortex method to the 
known value of downwash at the trailing edge. At a 
subsonic trailing edge satisfying the Kutta condition, 
— wjaU is unity; whereas at a supersonic trailing edge, 
— w[aU can be computed by the method presented in 
reference 1. 

Linearized theory neglects the effect on the downwash of 
the friction wake and the displacement and the distortion 
of the trailing vortex sheet. An experimental program, 
such as that reported in reference 10 for wings in subsonic 
flight, is ultimately required to determine the necessary 
modifications to linearized theory that will result in good 
agreement between theory and practice. 

SUMMARY OF ANALYSIS AND APPLICATIONS 

The perturbation field induced by a line vortex in a super- 
sonic stream and the downwash field behind a supersonic 
lifting surface have been examined to establish approximate 
methods for the calculation of supersonic downwash. 

An infinite line vortex of constant strength and slope, 
supersonically inclined to the free stream, induces no per- 
turbation field. A subsonically inclined line vortex has 
properties similar to those of a vortex in an incompressible- 
flow field. Bends- in a fine vortex induce infinite vertical 
perturbation velocities on the surface of the downstream 
Mack cone from the bend (except in the 2=0 plane). 

The downwash field behind a supersonic lifting surface 
differs from that behind a subsonic wing in several respects. 
For a supersonic lifting surface, discontinuities in shed vortic- 
ity occur at those points along a supersonic trailing edge 
where the plan-form slope is discontinuous and the local 
pressure coefficient is not zero. These discontinuities lead to 
singularities in the downwash field in the 2=0 plane. Also, 
the vertical perturbation velocities are logarithmically infi- 
nite on the downstream Mach cone from a wing tip formed 
by the intersection of a subsonic leading plan-form edge and a 
supersonic trailing plan-form edge. 

A bent-Iif ting-line method has been proposed for the solu- 
tion of the downwash field behind swept wings. When ap- 
plied to a triangular wing, this method gave results that were 
in very good agreement with the exact linearized solution for 
points near the line y= 0, 2=0 except for points within % 
chord of the wing trailing edge. 

An unbent lifting line (horseshoe-vortex system) has been 
proposed for unswept wings. This method was applied to 
determine the downwash behind rectangular wings with 
aspect ratios 6/e, of 2 and 4. Excellent- agreement was ob- 
tained for both aspect ratios by placing the lifting line at 
the X-chord point. 

Lewis Flight Propulsion Laboratory, 

National Advisory Committee For Aeronautics, 
Cleveland, Ohio, April 29 , 1949 . 





APPENDIX A 

SYMBOLS 


The following symbols are used in this report: 

) 

Cartesian coordinate system 


V, Vt 

2, 2, 

u 

v 

w 

i 

i 

k 

Z 

V 

r 

r 

& 

c . 
c p 

c , 

Cr 

E 


jperturbation velocities 
limit vectors 


( components of vortex vector (three-dimensional 
field) 

| components of vortex vector (vortex sheet) . 

wing span 

line integration 

local wing pressure coefficient 

effective chord of bent lifting line’ 

root chord of wing 

complete elliptic integral of second kind with 


modulus 


VK0 


Gi(m)= 
GM = 

m 

M 

m 

Vi = 

S 

u 


r fPr miXiY+tt-pmflZti 

-1 XiYt 


r< (Y t *+Z t *) 

spanwise lift distribution 

free-stream Mach number 

slope of line vortex or lifting line 


VA7-/W-0 2 Z t 2 ~ 

surface integration “ 

free-stream velocity (taken in x-direction) 
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V 

volume integration 


X ( = 

X — Xi 


Y< =- 

- y-Vi 


= 

Z—Zi 


a 

angle of attack 


p 

cotangent of Mach angle, Vil/ 2 — 1 


T 

spanwise distribution of wing circulation 


r m 

wing circulation at midspan 


6 

integration interval 


K 

circulation 


P 

free-stream density 


4> 

. perturbation-velocity potential 


A <t> = 

<t>r — 


a 

cross-sectional area of vortex tube 


O) 

resultant vorticity, (£ 2 -f ij 2 +f 2 ) w 
vfinite part of divergent integral 


§ 

line integral along closed curve 


Subscripts: 


a, b 

points of intersection of forward Mach cone 

with 


line vortex, lifting line, or edge of vortex sheet 

B 

bottom surface of z= 0 plane 


cp 

section center of pressure 


i 

1, 2, 3, ... n points on vortex lines 


l 

plan-form leading edge 


V 

variable of integration 


T 

..top surface of z=0 plane 


t 

plan-form trailing edge 


Superscripts: 



-value of . function at point of discontinuity 

when 


approaching from negative y-direetion 


+ 

value of function at point of discontinuity 
approaching from positive y-dircction 

when 



APPENDIX B 

UPWASH INDUCED BY LINE VORTEX 


The upwash induced by the line vortex segment of figure 1 is 


jS 2 k rn 

2r Jo [(1- 


m\y—mx) dy„ 


The upper limit is at the intersection of the line vortex 
(y„=mx*) with the trace of the forward Mach cone in the 
®=0 plane [( 2 — r„) 2 — )S 2 (y— y 0 ) 2 — jS 2 z 2 =0] and is therefore the 
appropriate root of the expression appearing in the denomi- 
nator of the integrand. The roots are 


P 2 m T }y 0 i +2m(p 2 my-x)y 0 -{-m\x 2 -p 2 y 2 -p 2 z x )] 3/i 

Vi) 


( 11 ) 


— m(p i my-x)±pm^(y—mx) 2 +(l—p 2 m 1 )z 2 


so that equation (11) may be rewritten 
-2x(l— j3 2 my /z 


dy. 


(B2) 


P 2 m i n(y-mx) W ~ ^ ( Jo l(y*— y»)(y»— Vo)] 3li 

From equations (6) and (7), 

| 7 =— J(0l)— <7. . (B3) 


where 



and 


m=- 

(y. 

C= lira 
v*-+y* 

2 A(y b ) 

jyt—y<> 


— Iim 1 

r • 

2 


2(y g +y»— 2y.) 


— y»)_ 


= 0 
Thus 


2 (ya+yt) 


, CB4) 

(3/a— ytf\y*yb 

Substituting equation (B4) in equation (B2) and solving for 
w yields 

it (y—mx)(p 2 my—x) 


w- 


2 ir (x 2 —p 2 y 2 ~p 2 z 2 ) l,2 [(y—mx) 2 -\-(l—p 2 m. 2 )z x l 
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(B5) 



APPENDIX C 

LOADING IN VICINITY OF WING TIP 


The nature of the loading in the vicinity of a wing tip 
formed by the' intersection of a subsonic leading edge with 
a supersonic trailing edge (fig. 6) is to be determined. 

According to linearized theory, u T in the region of a sub- 
sonic leading edge is singular of the form 


Substituting equation (C3) into the denominator of equa- 
tion (C2) yields 


n d(x t —x,) 
dr “ dy„ 

d v° Vy.[/(y»)] 


(C4) 


1 

•», r =-- -== - 
■S'X„—Xi 

where x t is the equation for the leading edge as a function of 
y„. This relation may be deduced from equation (11) of 
reference 9. The corresponding wing circulation is 

r=A <^=2 f — =4 A /a : t —Xi (Cl) 

JXl -yjXo—Xi 

The derivative of equation (Cl) is 

2 d(x,—xi) 

r (C2) 

“Vo -yjxr^xi 

where (x,—Xi) is the wing chord as a function of y„. How- 
ever, (x t —x t ) must be of the form 

Xt—Xi =yo\J(yo)] (C3) 

(where Lf(l/o)] y .=o^0, 00 ) ™ order to satisfy the restrictions 
that at the wing tip 

(a) the chord is zero 

(xt Xi) v °_. 0=0 

(b) the slopes of the leading and trailing edges are neither 

equal 



nor in the free-stream direction 



Equation (C4) indicates that dV/dy 0 is singular of order M at 
y 0 = 0 for the wing tip of figure 6. 
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